Two-Dimensional Weakly Interacting Bose Gas in the Fluctuation Region 
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We study the crossover between the mean-field and critical behavior of the two-dimensional Bose 
gas throughout the fluctuation region of the Berezinskii-Kosterlitz-Thouless phase transition point. 
We argue that this crossover is described by universal (for all weakly interacting \ip\ 4 models ) rela- 
tions between thermodynamic parameters of the system, including superfluid and quasi-condensate 
densities. We establish these relations with high-precision Monte Carlo simulations of the classical 
\tp\ 4 model on a lattice, and check their asymptotic forms against analytic expressions derived on 
the basis of the mean-field theory. 
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I. INTRODUCTION 

The mean-field (MF) approach to the weakly interact- 
ing Bose gas (BG) is a well established theoretical tool 
jl], H However, it is not adequate in the fluctuation 
region of the superfluid phase transition. The situation is 
most dramatic in the two-dimensional (2D) case, where 
the size of the fluctuation region, AT, is almost insensi- 
tive to the smallness of interaction [h = 1): 



AT/T C ~ l/\n(l/mU) 



(1) 



Here T c is the critical temperature, m is the particle mass, 
and U is the effective long-wavelength interaction con- 
stant. The regime of weak interaction corresponds to a 
small dimensionless parameter 



mU < I 



(2) 



which close to the transition point is equivalent to the 
condition nil <C T (n is the particle density). 

For the purposes of the present paper the microscopic 
origin of the effective interaction U is not important. 
However, to make connection between our results and 
realistic experimental systems we briefly review how U 
relates to the interatomic interaction potential V(r). The 
value of U corresponds to the pair vertex (the sum of lad- 
der diagrams) with typical external momenta ~ n 1 ! 2 . It 
can be written in a generic form as 



U 



V Q 



1 + (r7iVb/47r)ln(l/W 2 ) ' 



(3) 



where microscopic parameters Vq and d depend on the 
case. The simplest one is that of a weak short-ranged 
potential satisfying the condition V(r ~ ro) -C l/mr^, 
where ro is the potential radius. In this case, Vq = 
fV{r)drr, and with logarithmic accuracy assumed in 
Eq. (J3|) one may put d = tq provided nr^ <C 1 (when 
nr o ~ 1 the logarithmic term in the denominator can be 
neglected) . In the quasi-ZD system [|[ |(| , when the local- 
ization length of 3D atoms in the direction perpendicular 
to the 2D plane (axis z) is much larger than ro, one has 



first to average the pair interaction over a wavefunction in 
the z-direction, 4>o(z). Now, Vq = (Ana/m) f \4> (z)\ 4 dz, 
where a is the 3D scattering length. Introducing the 
localization length by l~ l — J \4> (z)\ 4 dz, we have Vq — 
Aira j 'ml z . With the same logarithmic accuracy, in Eq. 
d&l z for nil < 1 ( cf - Refs - i §)■ Finally, the case of 
strong short-ranged potential [which in the context of 
weakly interacting gas implies ln(l/mg) 3> 1] formally 
corresponds to the limit Vq — > oo in Eq. (||). In this case 
the effective interaction depends only on the parameter 
ln(I/ncP) with d m r (see Ref. 0). 

Because of weak log-dependence of the fluctuation re- 
gion on interaction one may wonder whether the MF the- 
ory makes sense at all in 2D (apart from the academic 
limit of exponentially small mU), and, if it does, then 
when. As a characteristic example of how problematic it 
is to reach the proper asymptotic limit, consider a dilute 
gas with very small nd 2 when mU ~ 4ir / \n(l / nd 2 ) <C I . 
With the same logarithmic accuracy it follows then that 
the critical point can be found without even resort- 
ing to the Berezinskii-Kosterlitz-Thouless physics || Q , 
T c « (27rn/m)ln _1 (f/TOf7). Meanwhile, a more accurate 
result for the critical point is (for future reference and 
convenience we write the answer for critical density as a 
function of temperature) [pi : 



n c = ^\n(^-) , e = 380±3. (4) 
27r V mU 



Obviously, an enormous value of £ makes it virtually im- 
possible to reach the limit of small U when £ can be 
ignored. 

There is, however, a very important point about the 
fluctuation region of a weakly interacting BG: In the limit 
of small U all l^l 4 models — quantum or classical, con- 
tinuous or discrete — allow a universal description J3|, [|. 
This observation follows from a simple fact that inter- 
actions are important only for long-wavelength compo- 
nents of the order parameter field, with momenta 
k < k c = m\J\JT <C fcr = V mT , and in this limit the 
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effective Hamiltonian is given by the \4>\ A model 

HM= /{^|VVf + |M 4 -MM 2 }dr, (5) 

where [i is the effective chemical potential. The micro- 
scopic physics of the model is important only at much 
higher momenta, k 3> k c , where the system behavior is 
ideal (in linear in U approximation) and thus may be 
easily accounted for analytically. 

This observation was successfully used in Refs. |[ [u]. [ll] 
(both in 3D and 2D) in the study of the critical point 
dependence on interaction. Same considerations apply, 
though, not only to the critical point itself, but to the 
whole fluctuation region around it, and one thus expects 
that, e.g., the superfluid density dependence on density, 
n s (n — n c ), or chemical potential, n a (fj, — (x c ), is also uni- 
versal close to the transition point and into the region 
where the MF theory takes over. 

The study of this universal behavior is the subject of 
this paper. We found that even for very weak interac- 
tion, say, mil ~ 0.01, the conventional MF theory result 
n s /n = 1 — T/T c may not be used since the fluctuation 
region is still of order T c itself. However, if T c is related 
to the density by Eq. (0), the modified version of the MF 
theory developed in this paper works remarkably well. In 
particular, the Equation of state and the quasicondensate 
density may be predicted very accurately up to T c (this 
is not true for the superfluid density, though). 

In Sec. we establish the universal form of the Equa- 
tion of state and the dependence of the superfluid density 
and quasicondensate density on chemical potential along 
with their asymptotic behavior away from the critical 
point. In Sec. ^ we describe the numeric model and 
the simulation procedure. Our results are presented and 
compared to the MF and Kosterlitz-Thouless theories in 
Sec. |y|. We conclude in Sec. [v] with discussing the ob- 
tained universal results in the context of quantum Bose 
gases. 



II. UNIVERSAL RELATIONS FOR WEAKLY 
INTERACTING V 4 MODELS 

The critical point of the BG is defined by Eq. (|J) and 
the corresponding relation for the chemical potential H 



mTU , 
Mc = In 



mil 



13.2 ±0.4 



which can be rewritten in the form 
jj, c = 2n c U H In 



mTU . 



(6) 



(7) 



Both n c and fj, c are model specific, and their values de- 
pend on the ultra-violet cutoff, , through the logarithm 
ln(fc*/fc c ); for the quantum gas fc* ~ hp- However, if the 
dominant MF type contribution to the chemical poten- 
tial, 2nU, is subtracted, the difference is an ultra-violet- 
cutoff-independent quantity. The same is true for the 



difference fi — \i c or n — n c . It seems natural then to 
introduce a dimensionless variable 



X = (fi — fi c )/mTU 



as a universal control parameter with the typical vari- 
ation across the fluctuation region of order unity. The 
Equation of state may then be written in the universal 
form as 



2nU - ^ 
mTU 



0(X), 



(9) 



where 9 is a dimensionless function. By subtracting crit- 
ical values from n and fi we can restate it as 

n-n c = mT\(X) , X(X) = [9(X) -8 +X]/2 , (10) 

with 



6(0) 



ln(£/U 



(11) 



From previous results |8) we have 6> = 1.07 ± 0.01, with 
the error bar being largely determined by the uncertainty 
in The X(X) function describes the so-called adsorp- 
tion isotherm which is relevant to the situation where the 
2D system is formed by atoms adsorbed on a surface. In 
the case of a trapped gas, the function X(X) describes 
the density profile of the gas in the hydrostatic regime 
(see the discussion in Sec. [V|). 

The behavior of the superfluid density is described by 
a dimensionless function /: 



n s = (2mT/Tr) f(X) 



(12) 



According to the Kosterlitz-Thouless theory |12| (see also 
below) , 



f(X^+0) — -> 1 + V2^X , 



(13) 



where k' is some constant to be defined numerically. 

Finally, in the superfluid region an important quantity 
is the quasicondensate density no, which we define by the 
relation 



M 4 



(14) 



The idea behind this definition is as follows. The notion 
of the quasicondensate || implies that the field ip has the 
following structure 



n e 



i<I>(r) 



(15) 
(16) 



where no is called the quasicondensate density, and ipi is 
the Gaussian field independent of ?/V Under these condi- 
tions, Q = Uq. [The Gaussian field ip\ obeys the Wick's 
theorem and thus does not contribute to Q.] Away from 
the superfluid region the notion of quasicondensate grad- 
ually becomes ill-defined, but the quantity Q is still of 
interest as a measure of local non-Gaussian correlations. 
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We will jargonically use the term 'quasicondensate den- 
sity' even well inside the fluctuation region, understand- 
ing by no the quantity \fQ. 

Since the MF theory result predicts no ~ n s , it is ap- 
propriate to characterize the dependence of Q on /i in 
close analogy to Eq. (|l2|): 



(17) 



Q = (2m7» g(X) 



The three functions — 6{X), f(X), and g(X) — com- 
pletely characterize system properties in the vicinity of 
the critical point 

We now turn to the MF and Kosterlitz-Thouless the- 
ories to establish asymptotic behavior of functions 6(X), 
f(X), and g(X). 

Asymptotic behavior at X — > oo. The notion of quasi- 
condensate is well-defined in this region and its density 
n a obeys a typical MF relation 



(n + 2n')U = n , 



where 



n 



(18) 



(19) 



is the non-quasicondensate part of the particle density. 
[One may worry how far the analogy between the gen- 
uine condensate and quasicondensate goes; the answer is 
that apart from the long-range order problem they are 
indistinguishable at the MF level, and our simulations 
confirm this assertion.] Comparing (JT^)- ( 19) with 
we see that 



n {X 



mT9{X) 



or 



g(X -> co) — (tt/2) 9{X) 



(20) 



(21) 



An explicit expression for the non-quasicondensate 
part is given by |l3|: 



d 2 k 



e(k) + n U - E(k) e(k) v E 



2E{k) 



E{k) 



(22) 



where e(fc) = k 2 /2m is the free particle dispersion law, 
E{k) — y/e(k)[e(k) + 2n a U] is the Bogoliubov quasipar- 
ticle spectrum, and ve — [exp[£'(fc)/T] — 1] _1 is the Bose 
distribution function. In the region of interest, 9 ~ 1, the 
first term in the integral is smaller than the second one 
by the gas parameter mU <C 1, and should be omitted 
fL6[. With the same accuracy, the second term yields 



n « -(mT/27r)ln(2n [//T) 



(23) 



With the help of Eq. ([2 
may be now written as 



the total density n = n! 



n 



26- -ln(20) 



mT 



1 



29 - - ln(2£0) 

7T 



(24) 



Substituting this relation into ( |To| ) we find the asymp- 
totic behavior of 9(X): 



9 -it- 1 In 9 -> X + 7r" 1 ln(2^) at X 



oo 



(25) 



To find f(X — » oo), we consider the standard expres- 
sion for the normal component density p| 



1 

2m 



d 2 k 
(2^P 



dE 



k 



(26) 



and compare it to the expression (|22|) for the non- 
quasicondensate part of the particle density. After in- 
tegration by parts and a straightforward algebra we find 
that up to higher-order in mU terms, n n — n' = mT/2n, 
which means that 



no — mT/2ir 



(27) 



and, accordingly, 



j) — » g(X) - 1/4 — » (tt/2) 9{X) - 1/4 . (28) 



It is important to note that while for obtaining asymp- 
totic relations (|2^) and (p7|)-(p8|) we employed the theory 
of the weakly interacting quantum Bose gas, the final re- 
sults are valid for any weakly interacting 2D system of 
the l^l 4 universality class, since these pertain to the uni- 
versal long- wave behavior of the system. 
Asymptotic behavior at X — > — oo. In the region X < 
the /-function is identically zero; the quasicondensate 
density is of no special interest in the normal phase 
[g(X) — * in this limit]. Hence, the only quantity we 
have to look at is the Equation of state, 9{X). Once 
again, we resort to the MF equation for the effective 
chemical potential fx' = ji — 2nU = —9mUT, and cal- 
culate the total density from the integral: 



d 2 k 



[exp(e/T + 6mU)-l] 1 



(29) 



-(mT/2n) ln(9mU) = n c 



mT 
~2T 



ln(0fl . (30) 



This expression may be immediately related back to the 
Equation of state (^|) and leads to the relation 



7T _1 In 9 



\X\ — 7T -1 ln^ at X — > — oo . (31) 



One has to understand the limit \X\ oo in the fol- 
lowing sense: it describes the system behavior close to the 
transition point but outside the fluctuation region. For 
the quantum gas (mU <C 1) it means 1 -C \X\ <C 1/mU. 
Of course, one may easily calculate system properties 
for any \X\ 1 using MF theory presented above, 
and, take care of the phonon contribution to the non- 
quasicondensate and superfluid densities at n ~ T/U 
[contained in the integrals of Eqs. ( p2| ) and (p6|)], or, in- 
stead of Eq. (|3(]), consider a more accurate expression for 
the dilute density limit n m -(mT/2n) In [l - e ~ emU ] to 
include the Boltzmann gas into the picture. Such obvious 
generalizations are not considered in this paper. 
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The vicinity of the Kosterlitz-Thouless point. The ther- 
modynamic limit close to the point of the Berezinskii- 
Kosterlitz-Thouless transition requires simulations of ex- 
ponentially large systems, and thus can hardly be treated 
numerically without renormalization group (RG) analy- 
sis of finite size corrections. Fortunately one can take 
advantage of the Kosterlitz-Thouless RG equations that 
describe the flow of the superfluid density n s (L) with in- 
creasing the system size L. In terms of dimensionless 
function ft — (n/2mT) n s (L) these equations are |12| 



6(X) 



dlni 



dy 
dlni 



= -y 2 l 



L ■ 



2(1 -h)y, 



(32) 



(33) 



where y(L) is the vortex-pair fugacity. By excluding vari- 
able y and integrating the remaining RG equation, one 
obtains the following relation 



F(f L2 ,f Ll ,K) = 41n(L a /Li) 
where F is defined as an integral 

dt 



F(a,b,K) = f 

J a 



t 2 (\nt- k) +t 



(34) 



(35) 



and k(X) is a size-independent (at k c L 3> 1) parameter. 

By performing large-scale simulations of systems with 
different sizes L\ < L 2 < L 3 < one may solve Eq. (|34|) 
for parameter k(X), verify that it is system size inde- 
pendent, and then determine the thermodynamic value 
f{X) = fh=oo{X) from the relation 



1// + 1h/ = k 



(36) 



that immediately follows from (p4|)-(p5|) at L — > oo. Since 
Eq. ( (36| ) has a root only at k > 1, we conclude that 
k = 1 corresponds to the critical point. In contrast to 
the superfluid density, k(X) has no singularities at the 
critical point X = and may be expanded into Taylor 
series 



k(X) ps 1 + k'X + 



(37) 



The solution of Eq. (|3^) for small X is then given by 
formula (13). 



III. NUMERICAL PROCEDURE 

Although all derivations presented in the previous sec- 
tion were done for the quantum BG, we expect them to 
be universal and apply for any model with effective long 
wave-length Hamiltonian (||) with small mil. Classical 
lattice algorithms are much more efficient than quantum 
ones and allow high-accuracy simulations of very large 
system sizes. Also, simulations of the classical lattice 




4 X 



FIG. 1: The Equation of state 8(X) compared to its asymp- 
totic large \X\ expressions, see Eqs. (G3) and (|3l[) 



model directly test the idea of universality, since they 
have to agree with all Eqs. (0), @, @, @, @. 

Our simulations were done for the simple square lattice 
Hamiltonian 



U 



(38) 



keBZ 



where ipk is the Fourier transform of the lattice field 
ipi, and E(k) = [2 — cos^a) — cos(k y a))/ma 2 is the 
tight-binding dispersion law; momentum k being defined 
within the first Brillouin zone (BZ). We employed the re- 
cently developed Worm algorithm for classical statistical 
systems |nj which has direct Monte Carlo estimators for 
all quantities of interest here and does not suffer from 
critical slowing down. 

We performed simulations for system sizes L = 64, 
128, 256, 512 and two values of interaction strength U = 
1/4 and U — 1/16 to eliminate finite-size and finite-C/ 
corrections to the results. Each quantity for each point 
in X was calculated with relative accuracy better then 
10~ 3 (down to 10~ 4 for smaller system sizes). Our final 
results for 6(X), f(X), and g(X) presented in Figs. ||, 
|3, and below are thus obtained with accuracy better 
than 1 % (the largest error bars of order 1 % are in the 
vicinity of the critical point where finite-size corrections 
are the largest; the errorbars for large \X\ shrink down 
to 0.3 -j- 0.1 %). Error bars are shown in all plots but 
typically they are much smaller than the point size. All 
the relevant data are mentioned in Table |. In Figures 
|l[|| the MC data are presented by dots and compared 
with the asymptotic analytic solutions shown by lines. 



IV. SIMULATION RESULTS 



5 



4 
3 
2 
1 

-1 



K(X) 



X 




FIG. 2: The X(X) function compared to its asymptotic large 
|X| expressions 
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FIG. 3: The quasicondensate density dependence on X com- 
pared to the asymptotic large \X\ behavior predicted by 
Eq. (pll) 



In Figs. [I], ^, |3|, and || we show our final results for the 
scaling functions. To compare them with analytic pre- 
dictions we first determine the value of the 9q parameter 
from the 9{x) plot to verify that it agrees with the result 
predicted by Eq. (pr|). We find that 



9 = 1.068 ±0.01 



(39) 



thus confirming the universality of parameter 9o- Know- 
ing £ M is all we need to handle the asymptotic behavior at 
large \X\. The best fits correspond to £ M = 13.4, which, 
within the error bars, coincides with Eq. (||). The agree- 
ment between the data and asymptotic laws Eqs. j25| ) 
and ( fjl] ) in Figs, [l] and || is remarkable. 

The same is also true for the quasicondensate density 
(see Fig. ||). The data perfectly agree with the idea of 
quasicondensate which is not entirely obvious at a first 



FIG. 4: RG parameter k(X) obtained from finite-size scaling 
of the data according to Eqs. ( pi) ) and (^5|) and fitted using 
linear in X expansion with dn/dX = 0.61 



f(X) 



1/f + ln(f) = l + K'X 
?0(X)-l/4 ,/* 



X 



FIG. 5: The superfluid density dependence on X compared 
to the asymptotic large \X\ behavior according to Eq. (pi[), 
and the Kosterlitz-Thouless theory for X < 0.25 according to 
Eq. © 



glance for the system without long-range order. We thus 
confirm that no, rigorously defined through the correla- 
tion function Eq. ( |l4| ) , plays the same role as the genuine 
condensate density in the 3D theory. 

As explained above, the analysis of the superfluid 
density data near the critical point is based on the 
Kosterlitz-Thouless RG equations. We found that for 
X 6 (—0.1, 0.25) the data scale according to the RG 
Eqs. ( J34| ) and (pq) with negligible finite-C/ corrections, 
and k(X) dependence is well described by the linear in 
X expansion Eq. (|37|), see Fig. ^. Having determined the 
first derivative of k at the critical point as 



dX 



= 0.61 ±0.01 , 



(40) 
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FIG. 6: Temperature dependence of n s /n for small mU = 
10"\ 10~ 2 , 1(T 3 , 1CT 5 . Points are connected with lines to 
guide the eye. The dashed lines are the mean-field theory 
results with exact relation between T c and particle density 



FIG. 8: Temperature dependence of no/n for mU = 0.2 de- 
rived from universal relations (filled circles) and simulated for 
the quantum lattice BG (open circles). Points are connected 
with lines to guide the eye 
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FIG. 7: Temperature dependence of no/n for small mU = 
10"\ 10~ 2 , 10~ 3 , 10" 5 . Points are connected with lines to 
guide the eye. The dashed lines are the mean-field theory 
results 



we proceed with the solution of the thermodynamic 
value of f(X) using Eq. (|36|). The results are plot- 
ted in Fig. [| along with the asymptotic large X law 
f(X) -> (n/2)6(X) - 1/4. The two limits match al- 
most perfectly around X = 0.5 and describe all the data 
points with exceptional accuracy. 



CONCLUSIONS 



The most striking result of this study is that the Equa- 
tion of state, 0{X), and the quasicondensate density, 
g(X), are predicted by the MF theory with accuracy of 



few percent all the way to the critical point. Even at the 
critical point the difference between the data and asymp- 
totic relations is barely visible. This outcome is quite 
unexpected because the superfiuid density does show de- 
viations from the MF theory for X < 0.5. 

Hence, from the side of the superfiuid phase, the 
boundary of the fluctuation region corresponds to X ~ 
0.5. To express this estimate in terms of temperature 
(density), we relate X and T at a constant density (X 
and n at a constant temperature). Equation ( |Io|) with 
n c from Eq. (Q) may be written as 

mT 2ir 

~ = ]n(£/mU) + 2ir\(X) ' 



or 



T c (n) 



n c (T) 



= 1 



2tt 



\n(£/mU) 



X(X), 



(41) 



and maps the control parameter X onto T/T c (n) or 
n/n c (T). This relation is not universal because T c , n c , 
and the r.h.s. of Eq. ( pl| ) depend on the ultra-violet cut- 
off. For X — 0.5 and, say, m[/ = 0.1 we find that the fluc- 
tuation region is roughly (T c — T)/T c ~ (n—n c )/n c ~ 0.3. 

It is instructive to compare results for the superfiuid 
density of the weakly interacting BG in the more conven- 
tional n s /n vs T/T c plot at a constant particle density 
which may be immediately obtained from the universal 
relations. The superfiuid density is given by 

±-™, mm J®%, m , (42, 
n nn m(£/mU) 

which along with Eq. ( ^l|) defines a parametric depen- 
dence of n s /n on T/T c . One may construct then a mod- 
ified MF theory result by substituting f(X) in this re- 
lation with n8(X)/2 - 1/4, see Eq. <M), and 0{X) from 
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Eq. (^5|). An analogous expression for no is obtained by 
replacing f(X) with g(X). 

In Fig. H we show the comparison between the MF so- 
lution and the data for small effective interaction mU — 
1CT 1 , 10~ 2 , 10~ 3 , 10~ 5 . We see that the fluctuation 
region is still of order T c even for mU = 0.01 — clearly 
that small coupling parameter may not be obtained by 
going to the dilute limit for the gas of hard-core particles 
and thus any realistic discussion of the experimental data 
should involve the proper description of the fluctuation 
region. Even for mU = 0.01 the familiar MF formula 
n s /n 1 — T/T c does not work at all — the slope keeps 
changing with T and the left-most slope does not point 
to n s /n = 1. However, if one uses an exact relation be- 
tween T c and n, then the MF approach described above 
is capable of reproducing the data for X > 0.5. 

The other remark concerns the asymmetry of the fluc- 
tuation region at T < T c and T > T c . The minimum in 
the Fig. [I] plot for 9{X) and the agreement with the MF 
laws already suggest that the fluctuation region is much 
broader on the normal side. When MF results for n s {T) 
in Fig. ^| are extrapolated to higher temperatures, the 
intersection with the temperature axis is still at sa 1.5 T c 
for mU as small as ~ 10 -5 ! We draw the same conclu- 
sion from the no/n plot shown in Fig. ^: The decay of 
the quasicondensate density uq = \/Q, which is a mea- 
sure of local non-Gaussian correlations, is quite extended 
into the normal state. 

Our universal relations are obtained in the limit mU <C 
1. From a practical point of view it is important to es- 
timate a typical value of mU at which higher-order cor- 
rections to our results become unimportant. To this end 
we note that if we plot the quasicondensate density as a 
function of T/T c for various values of mU, see Fig. |^ we 
will find that for mU = 0.1 the ratio no/n exceeds unity 
already at T k 0.4 T c . This unphysical result tells us 
that non-universal corrections for the quantum BG are 
not negligible even for mU ~ 0.1 (at least for the qua- 
sicondensate density; for numerical reasons they might 
be smaller for other quantities). We thus expect that 
universal expressions established in this study are likely 
to work without limitations only for mU significantly 
smaller than 0.1. To understand the situation with the 
quantum corrections quantitatively, we compare in Fig. 
g our results for no/n at mU = 0.2 with the previously 
reported |l5| results for the quantum lattice model [ fL8[ . 
The comparison suggests that at T ~ T c the quantum 
correction to the quasicondensate density is ~ mU (in 
relative units). The sign of the correction is negative, 
that is we are dealing with a quantum depletion of the 
quasicondensate. 

For helium films on various substrates jl9j [20|, and 
spin-polarized atomic hydrogen on helium film ]2lfl , the 
value of mU is most probably of order unity j^]. 

In the recently created quasi-2D gas of sodium atoms 
p2| , mV is of order 10~ 2 , and this system is supposed 
to be described by our results very precisely. In exper- 
iments with trapped gases the quantity directly relevant 



TABLE I: Final results (after taking care of finite-size and 
finite-(7 corrections) for the scaling functions Q(X), g(X), and 

fix). 



X 


0(X) 


o(X) 


f(X) 


-4.0056 


2 8363(3) 


2657(3) 




-3.0056 


1.9603(5) 


0.3094(4) 




-2.0056 


1.1472(6) 


0.3866(6) 




-1.5056 


0.791(1) 


4561(81 




-1.0056 


0.514(1) 


0.581(1) 




-0.7556 


0.434(3) 


688(31 




-0.5056 


0.442(3) 


0.869(2) 




-0.2556 


0.630(4) 


1.214(4) 




-0.1056 


0.885(9) 


1.560(7) 




-0.0556 


0.973(9) 


1.680(7) 




-0.0156 


1.041(9) 


1.774(8) 




-0.0056 


1.061(4) 


1.800(3) 




0.0044 


1.075(5) 


1.821(4) 


1.077(5) 


0.0444 


1.137(5) 


1.908(4) 


1.274(4) 


0.0944 


1.208(5) 


2.013(4) 


1.433(5) 


0.2444 


1.415(6) 


2.319(6) 


1.823(6) 


0.4944 


1.734(6) 


2.800(7) 


2.37(1) 


0.9944 


2.334(9) 


3.721(8) 


3.348(6) 


1.9944 


3.469(9) 


5.47(1) 


5.135(10) 


2.9944 


4.554(9) 


7.19(1) 


6.87(1) 


3.9944 


5.631(8) 


8.870(6) 


8.58(1) 



to the experimental setup is X(X) since it describes, ac- 
cording to Eq. (|lo|), the density profile in the trapping 
potential which is smooth enough to guarantee the hy- 
drostatic regime. In this regime the density variation 
over the mode-coupling radius r c ~ l/k c [the data of 
Ref. H suggest r c w 2/m(UT) 1 / 2 ] is small, and the coor- 
dinate dependence of density reduces to n = n(T, /x(r)), 
where n(T, u) is the homogeneous equation of state, 
u(r) = V^ x t(r) + const, and V oy ± is the trapping potential. 
It follows from Fig. [?] that for mU ~ 10~ 2 the size of the 
fluctuation region on the normal side is of order unity. 
Thus when the density at the trap center is tuned to the 
critical point, practically the whole density profile finds 
itself in the fluctuation region where MF equations do 
not work. 
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